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Abstract 

The prototype problem of a pendulum connected to a heat 
bath with strong dissipation is common for many physical 
applications. The corresponding equations of motion 
describe the current through the Josephson junction, and 
phase drift in a synchronizing auto-generator, is related to 
the problems of "stochastic ratchets", to the motion of 
Brownian particles in a nonlinear system with anisotropic 
spatially periodic potential. The thermodynamic flux in such 
systems is related to the (random) particle escape time. In the 
present work the first-passage time distribution and its 
thermodynamic interpretation has been used to describe 
these phenomena. In particular, according to the example of 
auto-generator oscillating system, it is shown that the 
unperturbed average first-passage time is related to the value 
of the thermodynamic kinetic coefficient. 
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Introduction 

The first-passage time, which is exit out of the system 
of quasi-equilibrium-state under the influence of 
fluctuations (Stratonovich, 1963, 1967), is important in 
many physical phenomena. This value is necessary to 
know in the theory of phase transitions, by the study 
of the dynamics of complex molecules, the chemical 
reaction rate, in the calculation of the coefficient of 
surface diffusion in semiconductors, as well as in the 
analysis of failure of radio and optical tracking systems, 
automatic tracking. In (Pontryagin and Andronov, 1933) 
Pontryagin equations were obtained for the probability 
of achievement, which is a classic work of Kramers 
(1940). Various aspects of this issue have been 
considered in (Maier and Stein, 1993). 

In (Stratonovich, 1995, Stratonovich and Chichigina, 
1996, Chichigina and Netrebko, 2001, Chichigina, 1999, 
2002) it was shown that the probability of the lifetime 


can be described by exponential distribution and the 
expressions for the average lifetime are obtained. In 
(Stratonovich, 1995, Stratonovich and Chichigina, 1996, 
Chichigina and Netrebko, 2001, Chichigina, 1999, 2002) 
a quasi-equilibrium Boltzmann distribution has been 
introduced (Chichigina and Netrebko, 2001, 

Chichigina, 1999, 2002) and it is shown that the 
entropy decreases in the measurement of the lifetime 
(Chichigina and Netrebko, 2001, Chichigina, 1999, 
2002). These results can be used to support the 
approach of (Ryazanov, 2004, 2006, 2011), in which the 
time to achievement of level is considered as 
thermodynamic variable. In contrast to (Stratonovich, 
1963, 1967), a lifetime without impacts is separated, 
which is part of the generalized structure factor, and 
the value containing the impact on the system and 
thermodynamically conjugate to thermodynamic 
quantities of time to achieve level. The method of 
nonequilibrium statistical operator in (Ryazanov, 2001) 
is interpreted as an average of quasi-equilibrium 
statistical operator on the distribution of the lifetime of 
the system. In (Chichigina and Netrebko, 2001, 
Chichigina, 1999, 2002) a quasi-stationary distribution 
is obtained. Distribution obtained in (Ryazanov, 2004, 
2006, 2011), after the integration over of lifetime is 
stationary and the possibility of exit out the system is 
taken into consideration. In contrast to the approach of 
(Chichigina and Netrebko, 2001, Chichigina, 1999, 2002) 
in (Ryazanov, 2004, 2006, 2011), the smallness in 
lifetime rather than the value, as well as the 
insignificant effects on random variable of lifetime is 
proposed. Therefore the condition that the average 
lifetime must be much larger correlation times, 
mandatory in (Chichigina and Netrebko, 2001, 
Chichigina, 1999, 2002), not necessarily in the approach 
proposed in (Ryazanov, 2004, 2006, 2011). In this paper 
the results of (Ryazanov, 2004, 2006, 2011) in the fourth 
section shall be applied to physical problems described 
in the second section. In the third section, a 
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thermodynamic interpretation of lifetime is provided. 


Phase Synchronization 

In this section, the synchronization of the generator in 
the presence of noise is in consideration (Stratonovich, 
1963, 1967, Tikhonov and Mironov, 1977). Effective 
means to achieve a more stable oscillation frequency is 
the external clock, that is, to supply voltage for the 
clock generator from another generator, which has a 
greater degree of frequency stability, while with a low 
power. An entrainment of the oscillation frequency 
and its stabilization are in it, however, the presence of 
fluctuations leads to a small deviation from the phase 
of its synchronous value, and to a large deviation of 
the phase, comparable to n/2, as well as the phase 
jumps by an integer number of periods. These jumps 
failing to be compensated, are accumulated and lead to 
the fact that the average frequency of oscillation does 
not coincide with the clock frequency, as well as to the 
fact that there is an irreversible phase diffusion, which 
gives the diffusion frequency instability. 

If the intensity of the noise is small 

K « ectkAo 1 , (1) 

where co c is frequency of synchronization, the small 
and the relative variation of the amplitude 

((A-a, 2 ))«A> 2 ' (2) 


the exact value of the amplitude A incoming in the 
simplified equation for random phase <p can be 
replaced by an approximate value of Ao and the 
equation for the phase of the form is obtained 


dcp 2A% 

— = A 0 -Asin <p — — 
dt A 


( 3 ) 


In (l)-(3) A 0 = (coo 2 - o>?)l(Oc * ox -coc is the value of the 
initial mismatch of generators co, and coo=(LC)~ m is 
resonant frequency of master clock, in which C is 
capacity, L is inductance, R is resistance of circuit of 
simplified diagrams of oscillator, M is electromotive 
force induced circuit due to the presence of inductive 
coupling, S is slope of cubic lamp characteristics, A is 
amplitude. Am is amplitude of sync harmonic signal 
Amcoscot ; while the formula Ao= 4(MS-RC)/Myi is 
stationary value of the amplitude, in which ratio yi 
describes the extent of the cubic dependence, 
s=coo 2 (MS-RC)/ coc, and K is correlation parameter of 
random forces £ (intensity of the noise) from ratio of 
<^^>=K§(t), A=(OcE/Ao which is band of 


synchronization (retention), E=coo 2 MCUo, Uo is the 
amplitude of the sync voltage. 

Equation (3) is applied if (1) and 

E/Ao « s, ecoc » Ac, (4) 


A c= cocEHAo, when linearization of the amplitude and is 
not allowed linearization of the phase. Equation (3) 
corresponds to the Fokker-Planck equation of the form 


dco(cp) _ S[(A 0 - Asin<p)<y(y>)] K d 2 co(cp) ^ 
dt dcp 2Ag 2 dcp 1 


where (A <p) is probability density of the distributions of 
the random phase values cp. Note that equation (3) and 
(5) arise in many physical problems. Thus, in 
(FeigeTman and Tsvelik, 1982) this equation (with 
other symbols and, respectively, with different physical 
meaning of the parameters) describes the motion 
placed in a thermal reservoir pendulum with strong 
dissipation under the action of gravity and torque. In 
(Ivanchenko and Zil’berman, 1969, Ambegaonar and 
Halperin, 1968) this equation describes the current 
through the Josephson junction with a large dissipation. 
In addition to maintenance phase of clock oscillator 
(Stratonovich, 1963, 1967) it describes the care phase of 
the laser (Haken, 1980). In recent years increasing 
attention of researchers has been paid to the so-called 
"stochastic ratchets" (Magnasco, 1993) under the 
understanding of which is the motion of Brownian 
particles in a nonlinear system with a spatially periodic 
anisotropic potential. In such systems, when exposed to 
different types of fluctuations, it may cause the 
probability flux along the spatial coordinates. In 
(Postnov et al. 1996) the mathematical model of the 
radio system phase-locked with an appropriate choice 
of parameters is reduced to one-dimensional model 
with a spatially periodic sawtooth potential. 
Experimental results in a one-dimensional 
approximation showed the effect of emergence of a 
probability flux. In this case two different mechanisms 
for its occurrence have been implemented. The 
experimental setup corresponding to a non-zero 
probability flux drift mutual phases of the two 
oscillators, the direction and magnitude are determined 
by the parameters having affect on the installation of 
signals. 

Stationary solution of (5) is written in (Stratonovich, 
1963, 1967, Tikhonov and Mironov, 1977) as 

<p+lK cc.\ 

co st (cp ) = N~' exp{D 0 p + Dcosp) J exp {-D 0 i// - Dcosi//}dip W 
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N=47r 2 exp{-xD 0 }/[iDo(D)/ 2 , where Ldo(D) is tabulated 
Bessel function with imaginary argument and 
imaginary index, Da=Am 2 /AK, Do/D=A 0 /A. The 
parameter D characterizes the signal to noise ratio in 
the band of synchronization, and D 0 is the value of the 
relative initial detuning. In (Tikhonov and Mironov, 
1977) the value 

, _ A ,(±>. _ Aexp{(±);rD 0 } (7) 

J (±) — 1 — i ,2 

4x 2 d\I iDo (D)\- 

is introduced. That is average number of jumps of the 
phase difference by 2 n per unit time in the direction of 
its increase (decrease). The value 
J=Gi=Ni + -Ni~=Ash( 7rDo)A7r 2 D fliDo(D)f=A(psh( xDo)/27r 2 Do/li 
00 (D )/ 2 characterizes the average flow, and the value 
Ni=Ni + +Ni-=cth(ftD 0 }Gi- the average diffusion. In 
(Stratonovich, 1963, 1967, Tikhonov and Mironov, 1977) 
it was defined as the times of the transition of phase 
from the vicinity of a stable equilibrium state in the 
neighboring states of stable equilibrium. So for the 
average time Ti of first passage, the nearest stable 
equilibrium states from the initial steady state in 
(Stratonovich, 1963, 1967) the expression are written 

_ 4,- I) exp { zl).. )'/ j, ( I) l sh(nD 0 ) 1 _ T _ 4n 2 D\l i0 (D)f , 

1_ A(l + exp{2ffD 0 }) “ ch(nD 0 )J ” A^’ °” ,|Dg=0 ” 2A 


wells is small flows themselves /+, J- can be large, and / 
« /+. Conversely, if the system is highly 

non-equilibrium, then conversely /- « /+ and / ~ /+. The 
flow Jo (or the time to reach To) is, therefore, an 
inherent characteristic of the system, a kind of 
potential to respond the force. It seems to be similar to 
the kinetic coefficient. You also need to find out what 
time scale is in consideration. The consideration with 
streaming necessarily requires that the time scale has 
been much more time taken by the particle in the well 
(only then can we even talk about of the flow). 
Otherwise, the flow doesn't exist as statistical value. 
Thus, the formula for the frequency hopping in 
equilibrium is obtained by setting Do = 0 in (7) 


T, = Ah = - 


4tt 2 D\I 0 (D)\- 


( 9 ) 


this "equilibrium" frequency of hopping. Power is 
indicated here as Do, and D is parameter such as 
temperature. Now the same expression (7) in the 
presence of power should be paid attention to find 
uncompensated flow J = Ni + -Nr. Note that for small 
forces Bessel function is almost the same as Io 
(compared to the contribution of the leading 
members). 


So, it can be written 


( 8 ) 

To is mean time to reach the boundary in the 
equilibrium state. 

A Thermodynamic Approach to the 
Description of the J umps of the Phase 
Difference. 

It should distinguish between the different meanings 
of "flow." Thus, there is a "probability flux" out of this 
potential well into the next, as it exists in the 
equilibrium and in a non-equilibrium state. It is also 
associated with the time of reach of the boundary. For 
example, in the equilibrium state it is (average) flow Jo 
and the mean time to reach the boundary To=2/Jo=2/No 
(Magnasco, 1993) (symmetrical with respect to the left 
and right outputs). In the non-equilibrium state (at 
nonzero pumping Do) it is frequency hopping (flows) 
to the left and to the right /+ (=Ni + ) and J- (=Nr). And 
then there will be an uncompensated flow /= /+- J- 
(Gi=Ni + -Nr). It is this flow (and not /+) that is a 
thermodynamic quantity. The relationship between the 
thermodynamic flow J and J+, J- can be dependent on 
the magnitude with different degree of impact. For 
example, if the difference between the depths of the 


J ~(exp { nDo } - exp 1 -ttDo}) Jo, (10) 

that is, the thermodynamic flux / is expressed through 
"internal characterization" Jo. In the spirit of linear 
thermodynamics the exponential must be extended in 
(10) and left linear member the total acquisition is that 

J ~X Jo = 2X / To, (11) 

X is the power (in this case 2 nDo), and the coefficient of 
proportionality is related to the achieved time. That is 
(equilibrium) first-passage time similar to the kinetic 
coefficient. Since the results are utilized (for 
synchronizing the generator or for a nonlinear 
oscillator), it can be said that it is confirmation of the 
theory (Ryazanov, 2004, 2006, 2011) on the microscopic 
level. 

Heuristically (without specifying the model) the 
following consideration can be cited the general 
formula for the internal production of entropy is 
expressed as 

8S=JX. 

As in conventional linear thermodynamics, X can be 
expressed from (11), and we have 

SS= J 2 To/2. (12) 
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The positivity rate of entropy provides a positive 
lifetime. Thus, the unperturbed mean time is to reach 
the boundary maps to the kinetic coefficient in the 
approximation of linear nonequilibrium 
thermodynamics. 


Distribution with a Lifetime 


In (Stratonovich, 1963, 1967) the distribution is 
obtained 


co(cp,r) = ®[(^)exp{-Ar}, A = (r') ' / 

(13) 

(where r is time to reach the absorbing boundary), the 
corresponding distribution of (14) obtained in 
(Ryazanov, 2004, 2006, 2011). The value of r in this case 
serves as a lifetime (Stratonovich, 1963, 1967), and A 
plays the role in the conjugate of the thermodynamic 
quantity. Distribution (13) in (Stratonovich, 1963, 1967) 
corresponds to the "stable decrease" of the number of 
remaining in stationary area of the phase space of the 
representative points. In (Stratonovich, 1963, 1967) the 
value A is expressed in terms of the intensity of jumps 
of phase difference (7) and the average lifetime, A = y= 
Jch(7iDo)lsh(7iDo), where / is average flow of the second 
and third Sections, and equating A=y is based on 
equating the expressions (13) and (14) a different 
expression for y is obtained that does not use these 
assumptions. 

In (Ryazanov, 2004, 2006, 2011) the distribution 
containing lifetime (instead of r in (13)) is introduced. 
The joint probability density for the energy values of E 
and the lifetime is 


P(E, T) = 


- pE ~ rT co{E, T) 


z{p, Y ) 

Z(J3,y ) = J e~ pE ~ yE dz = // dEdru)(E,r)e-l 1E -v r 


(14) 


where z is dynamic variables (z =qi,...,qN / pi,...,pN), oXE, 
E) is generalized structure factor, the number of phase 
points between E,E+rf£; r,r+dr equals oXE,F dEdr 
(Ryazanov, 2004, 2006, 2011, 2005, 2012). 

Integrating (14) by r , the stationary distribution for E 
of the form is acquired 

e -P E re 

p(E) = \P(E,T)dr=— — je-^C^Ddr. (15) 

J Z{P,y) J 0 

A similar result was obtained in (Chichigina and 
Netrebko, 2001, Chichigina, 1999, 2002). However, in 
(Chichigina and Netrebko, 2001, Chichigina, 1999, 
2002), the distribution of the form (13) is obtained. 


which is, in general, not stationary. In (Ryazanov, 2004, 
2006, 2011) the generalized structure factor gives 
meaning to the stationary joint probability density for 
the E and F We write 

®(£,r) = ®(£')® 1 (£', r ) = co{E)Y j R k f k {r,E), (16) 

k = 1 

In the expression (16) it is assumed that the system has 

n classes of ergodic states; in which is the 
probability of the system in fc-th class of ergodic states, 

and f k <T.E) is density of lifetime r distribution in 
this class of ergodic states. As a physical example of 
such a situation (characteristic of metals, glasses) the 
potential functions of many complex systems can be 
specified (Olemskoi, 1993), including the potential of 
the system discussed in Section 2. Potential minima 
correspond to metastable phases, disordered structures, 
etc. In such systems the phase space is divided into 
isolated regions, each of which corresponds to a 
metastable thermodynamic state, and the number of 
these regions increases exponentially with increase in 
the total number of particles (or quasiparticles) 
(Edwards, 1994). Quasithermodynamic theory of 
structural transformation of the alloy Pd-Ta-H, based 
on this model, has been developed in (Avdyukhina et 
al. 2002). Expression (16) is applicable to the 
description of such systems, the explicit form of the 

distribution fk is defined in (16), namely the gamma 
distribution 

f k (x) = } \ x ai ~'e~ xlbt , x>0, f t (x) = 0; x<0 

r k)V‘ 

J 0 ” e-Y« x f k (.x)dx = (1 + Ykb k r ak (17) 

((a) is gamma function, bk are parameters of 
distribution f k (F=x)). Substituting (17) in (15) - (16), we 
obtain 


e -fiz “ « 

_ 6)(E)\YRJ k (r,E)e-' tT dr = 


-PE „ 

——aXEj^Rpx+rA) 

2(P,y) F! 

From (14) and the fact that dlnZ/dy/p = - <F>: 

</>*- = a k b k ; r 0 K =<rp> Kjr a=a K b v 

1 + rA 

l + y K b k= ; 


(1+ykbk) ak =exp{-akln( T ot )}. 


(18) 


(19) 
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In many cases, instead of the gamma distribution (17) 
exponential distribution can be used with a*=2. When 
ak=l and b K = To K in (Ryazanov, 2005, 2012) the 
distribution is obtained 

p(E) = 

s =1 U = ggo?) 3=1 R* [1 + f k exMPPV) 0k r 

where Pk, Vt, Tk are pressure, volume and temperature, 
respectively, in the fc-th class of metastable states, and 
index 0 at the bottom indicates the absence of 
disturbances, when ytc=0, 7a=7, is intensity of the flow of 
energy in the system (subsystem), which is equal to the 
dynamic equilibrium of the output intensity. This 
distribution tends to Gibbs at y—>0. For large values of y, 
p(E)-tf). At ymX, 

e -fB 

p(E)~ „ a , m( ~E)yR, : c -*j5P 'totEWft, Vk=Vk/V, u=E/V. 

Z{P,y) tf 

As a physical example of the value y/X consider its 
value for the random generator synchronized phase 
(Section 2), where (Stratonovich, 1963, 1967, Tikhonov 
and Mironov, 1977), average time to reach (in the 
approximation of the initial steady state) is defined by 
(8), where the value of Do is related to disturbance y, 
<rp>=Ti (13), ro=To=l/2A, 2 = tv = — is the 

° 4/r 2 Z)|/ 0 (D)| 2 

intensity of the phase jumps at Do = 0 (9). From (19) at 
«=1 , R>=b,. 1 1 ..Then 

7 j_ [l + exp{2;rD 0 }]|/ 0 (D)| 2 

1 1 ( r r) ( r <>) exp{^ 0 j|/, Do (D)| 2 -2 

nDg [exp { nD a } + exp { -flD 0 } | ' 

sh(nD a ) {1 + 2D 0 2 / 0 - 2 (Z))£ (-1)" /„ (£>) /(n 2 + D 0 2 ) } - 2 

n = 1 

This quantity tends to 0 by Do—XJ and to co as Do by 
Do—>co (although Do takes limited value). As the value 
y (20) is known, the distribution of the form (14), (15) 
can be used to solve all sorts of non-equilibrium 
problems, for example, for the thermodynamic 
description of the jumps of the phase difference 
(Section 3). 

Conclusion 

Thus, the use of the lifetime as a thermodynamic 
variable allows us to obtain a stationary 

nonequilibrium distribution for "primary" variable 
energy (in our example - the random phase). These 
distributions can not only describe in detail the 
processes of phase synchronization, but also a number 


of other phenomena mentioned in the Introduction. 

Comments have been made on the time scale and the 
behavior of the average time to reach the borders has 
been qualitatively discussed. Relative to the average 
time of exit out of the potential well, there are two 
different time scales. Quasi-equilibrium in a potential 
well is not the same as the general equilibrium when 
the average is taken over the set of potential wells (or 
periods of employment in terms of queuing theory). 
The exit time is large for a number of systems with 
small fluctuations, and these two scale are the same. 
However, in small systems (or systems with large 
flows), this is not so. Ordinary thermodynamics, 
Boltzmann equation and thermodynamic equation of 
Maxwell describe great times together with the 
assumption that exit time is also very large. In the case 
where it is essential to record the value of periods of 
employment, it is not - here for all the time dynamics 
lying on a small time scale, and a quasi-equilibrium is 
established on time scales shorter than the output from 
the potential well, while the global equilibrium is 
established for much larger lifetimes. Therefore, in 
such situations, it is difficult to find evolution 
equations. The dependence on memory can be on a 
scale larger than the time to reach. As in (Stratonovich, 
1963, 1967, Tikhonov and Mironov, 1977, Stratonovich, 
1995, Stratonovich, Chichigina, 1996, Chichigina and 
Netrebko, 2001, Chichigina, 1999, 2002), one can 
imagine the process of hopping between potential 
wells as Brownian motion, but the memory may be 
housed in the complex structure of individual well. 
This dependence of the memory can be neglected in 
the simplest cases. Accounting random time spent in 
the potential well leads to a semi-Markov processes (if 
these times of stay are not exponentially distributed). 
The assumption of exponential distribution 
(Stratonovich, 1963, 1967, Stratonovich, 1995, 

Stratonovich, Chichigina, 1996, Chichigina and 
Netrebko, 2001, Chichigina, 1999, 2002) leads to the 
Markov nature of the process. In the general case the 
distribution for the first passage time can be written 
with the dependence on the initial position. These 
distributions will be more complex than the 
exponential distribution, and the result will be a 
semi-Markov process, that is, it already has a memory. 
The average time to reach is almost constantly in the 
potential well, but it changes abruptly near its 
boundaries. If the potential well is deep enough, the 
corrections are of higher order. However at high 
temperatures or at high thermodynamic forces, the 
amendment on the borders can be compared with the 
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behavior in the depths of the pit. 
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